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ABSTRACT 


A general theory for indicial potential- compressible 
aerodynamics around complex configurations is presented. 

The motion is assumed to consist of constant subsonic 
or supersonic speed for tfO (steady state) and of small 
perturbations around the steady state for t>0. Using 
the finite-element method to discretize the space problem 
one obtains a set of differential-difference equations in 
time relating the potential to its normal derivative on 
the surface of the body. The aerodynamics transfer function 
is then obtained by using standard method of operational 


calculus . 
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INTRODUCTIO N 

Considered in this note is the problem of unsteady 
subsonic and supersonic potential aerodynamics for an 
aircraft having arbitrary shape. The' motion of the 
aircraft is assumed to consist of small perturbations 
with respect to the constant speed motion. The objective 
of this paper is to describe the time functional relation- 
ship between aerodynamic potential and its normal derivative 
(normal -wash, / 3 ■/) ) in a form which can be easily used 

for the flight-dynamics analysis in both time and frequency- 
domain. The finite element method is used for space- 
discretization. 

The analysis presented here is based upon a new 

1/2 

integral formulation, derived by the author, which 
includes completely arbitrary motion. However, the 
numerical implementation (Refs. 3 and 4) is limited to 
steady and oscillatory flows. On the other hand, in order 
to perform a linear-system analysis of the aircraft it 
is convenient to use more general aerodynamic formulations 
i.e. fully transient response for time-domain analysis, 
and the aerodynamic transfer function (Laplace's transform 
of the fully unsteady operator) for frequency-domain 
analysis (see for instance Ref. 5). Consistently with 
this type of analysis the unsteady contribution is assumed 



2 


to start at time t=0, so that for time t<0 the flow 
is in steady state. Furthermore the motion of the 
aircraft is assumed to consist of small (infinitesimal) 
perturbation around the steady state motion. Since the 
initial work by Wagner on unsteady incompressible two- 
dimensional flow, several problems have been considered. 
Detailed analysis of the various methods available are 
given for instance in Refs. 5 and 7. Finite wings can be 
solved only for particular planform (such as elliptical 
wings) for subsonic flows, and for general planform in 
supersonic flows. Strip theory is used for slender wings 
in subsonic flow. 

For subsonic and supersonic flows around arbitrary 
complex configurations no tool has been available for 
either time- or frequency-domain analysis. Such a tool 
is presented in this note, for both time and frequency 
domain. For the sake of brevity, only the relationship 
potential-downwash is outlined since the full transfer 
function (generalized force per unit generalized coordinate) 
can be easily obtained by including the relation downwash- 
generalized coordinate (boundary conditions on the body) 
and the relation between potential and generalized forces. 
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SUBSONIC INDICIAL AERODYNAMIC S 

Consider the subsonic case first. Within the 

small perturbation assumption, the motion of the surface 

of the aircraft with respect to a frame of reference 

traveling at uniform subsonic speed with respect to the 

undisturbed air can be assumed to be negligible. Thus 

the Green theorem for the equation of the aerodynamic 

1, 2 

potential is given by 

An E(P<) £(?, T ) - O]® j-dl + 




dt 

R 3H 


(1) 


where Z. is a surface surrounding the body and the wake, 

\]/ c 3 /3N is the normal derivative (normal-wash) of i 

on £ , N is the normal to the surface V 


z -- ?/i ; t, « si a 


and 


( 2 ) 


R = ffx-xj 2 + (y-y ¥ ) z + (z -z x f] 


( 3 ) 


while 


[ l 


e 


J T. r T - 0 
1 


( 4 ) 


* 
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where 

0 = A1 (V-X*) +* (5) 

is the time necessary for a disturbance to propagate 
from P to P* . Furthermore 

EE. = 1 o i ci&. 2. 

- i/2 on 1 

Z 0 t r\ £ f <£e 21 C6 1 

In Eq. 1 the surface 2! is assumed to be fixed with respect 
to the frame of reference. However, the effect of the 
motion of the surface is retained in the boundary condition, 
which gives the normal-wash approximately as 



(where ^j(pxJxO is the equation of the surface 2E ). 

This approximation is consistent with the hypothesis of small 
perturbation which has been invoked also in writing Eq. (1). 
Also consistent with the hypothesis of small perturbation 
with respect to the constant-speed motion, is the assumption 
that the surface of the wake is the one of the steady-state 
case. 


It may be noted that, because of the linearity of the 
problem, the state-steady contribution can be subtracted 
from Eq. (1) . Therefore in the following it is understood 
that ^ and ^ are the unsteady parts of the potential and 
the downwash, which (in line with the concepts of operational 
calculus) are assumed to be identically equal to zero 
for T < O : 



- 5 - 


$ : O ; I'iO (t < O) (8) 

In order to understand the nature of the aerodynamic 
operator, Eq. (1), it is convenient to isolate the 
contribution of the wake. This yields 


4»E(P.)4>(P„,T) 



is the (open) surface of the wake, and is the 

potential-discontinuity across of the wake, evaluated 
in the direction of the normal (i . e. if the 

upper normal is used) . It should be noted that the value 
of is not an additional unknown, since it can be 

easily shown from the Bernoulli theorem that (see 
Appendix B) 

A$(?,T) = , T-n) (10) 

where H is the nondimensional time necessary for the 
vortex-point to travel (within the steady flow) from 
the point, P TE (origin of the vortex-line at the trailing 
edge) , to the point P. 

If small perturbation hypothesis can be used for the 
steady state flow, then TT can be approximated by 

71 - P* (x-x TE )/' n U1) 
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Equations (9) and (10) fully describe the problem of 
linearized unsteady subsonic potential aerodynamics 
around complex configurations. In order to solve this 
problem, it is necessary to obtain a numerical approximation 
for Eq, (9) . This can be obtained by using for instance 
the finite-element method. While other methods can be 
used, the finite element one appears to be at present 
time the only method sufficiently general and flexible to 
be used here. Consider the integrals on first. 

Using a typical finite-element representation, it is 
possible to write 

^0 /T -o;= £ -©) m l (p) 


<£ (p, t- q) = 2 © u ) m^(p) 


a2i 


where is the total number of nodes on the body, 

f;(T-6>).»a 4JT.0J are time dependent values of 

and 4* at the node, L, at the time T-0 (where c5 is 

w * l 

the disturbance-propagation time from P* to P ) ? 

furthermore M (P) and N (P) are prescribed global 
■U L 

shape-functions, obtained by standard assembly of the 
\ 

element shape-function (see for instance Ref. 8) . 

For instance for the hyperboloidal element of Ref. 4 the 
element shape functions are 


M/ -=n; -.!(> -i/qo ■)/■),) 


(13) 
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where £ =t i and i] - ± 1 define the locations of the 
corner, A , of the element, E. 

Next consider the integrals on the wake. In order 
to facilitate the use of Eq. (10) , it is convenient to 
divide the wake into strips defined by (steady-state) 
vortex-lines emanating from the nodes on the trailing 
edge. The strips are then divided into quadrilateral 
elements. The potential discontinuity can then be 
expressed as 

A ch (?, X- e), t * £( X- ©,) LJ-P) il4) 

where N is the number of nodes on the wake is the 

° 'IT 

value of at the N-th node on the wake at time 

T- © (where 0 is the propagation time from P to P ) 

N W * f 

and L^(P) is the global shape function relative to the 
N-th node of the wake. Note that according to Eq. (10.) 

where M=M(N) identifies the trailing-edge node which is 
on the same vortex-line as the node N. Furthermore H .. 

/V 

is the time necessary for the vortex-point to be convected 

from the trailing-edge node M to the wake-node N. It 

may be worth noting that & z. S - cb , where 

1 f\ 'U 

L u and identify the upper and lower trailing-edge nodes 
on the body corresponding to the M node on the 
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trailing edge. Therefore it is possible to write 

= Z $ cb 

W(fl } NL L ( 16 ) 

where r 5 -1 ( if L identifies the upper 

hi ^ A.' L- > 

(lower) node of the body corresponding to the N-th node 
of the wake , ar>A o otherwise . 

Combining Eqs. (9, 12, and 14) one obtains 

2E(P V ) 4>(IJ) = T \ t L (T- ©J + 



C 18 ) 
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and according to Eqs . - (15) and (16) 

t S -nj 


(19) 


Next consider, in particular, that P* coincides with 
the node J of the body. In this case E= 1/2 and, using 
Eq. (19) , Eq. (17) reduces to 

= ? V t ( T ' Q «.) + l C ]L 4>,&- &J) 4- 

? \ i t< T - ®-) * H 


Z. X G S — <fv f 1 “ - ^0 

jtj hi l ^ ^ ^ ^ / 


t. 


(2 0 ) 


where 


(b„ , C» , t>„ , t* , , S J ! fc , c , \ , r. ,s, , q) 


( 21 ) 


Equation (20) indicates the nature of the aerodynamic 
operator relating potential and normal-wash as obtained 
by using finite-element representation to discretize 
the spacial problem. The operator is a linear differ- 
ential-difference operator to which the methods of 
operational calculus can be applied. However before 
considering the Laplace’s transform of Eq. (20) it is 
convenient to make some remarks about the contribution 
of the wake. It may be noted that, according to Eq. (8), 
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<|> is identically equal to zero for T<0 therefore 
according to Eq, (15) 


A <3? , 3 0 

Al 




Hence if the analysis is limited to T < T 


max 


( 22 ) 


the 


contribution of the elements with H > T is identically 
equal to zero. Therefore the wake can be truncated to 
eliminate these elements. It may be noted that those ' 
elements would contribute to the transfer function and 
thus to the transform of $ but not to the final solution 
in the time domain for T < T . The advantage is not 
only that less computational time is used .(since less 
elements are required) but also that the problem of con- 
vergence connected with the infinite wake (factors e 
with Real(s) <0 and D oc ) are bypassed.* 

rJ 

Next by taking the Laplace transform of Eq. (20) 
and solving for one obtains 




(23) 


* A correct analysis implies the evaluation of the 
limit of the present analysis as the number of elements 
on the wake goes to infinity. 
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l 



Equation [ 22 ) indicates that the matrix [ Aj^-3 is the 
desired subsonic matrix transfer function relating the 
transformed vector of the potential £ to the trans- 
formed vector of the normal-wash 
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SUPERSONIC INDICIAL AERODYNAMICS 


In this section the formulation for the supersonic 
case is briefly outlined. For simplicity only supersonic 
trailing edges are consideres so that the contribution 
of the wake can be ignored.* Under small perturbation 
assumption the Green theorem for potential supersonic 
flow is given by 

4nE(PJ 4(t,T) * - 41 + 


-£ 

©’ 


R 




(25 ) 


where " , 5 r = d&/dM e ’ ( is the conormal derivative^) 

is the conormal-wash, 

X=^/(g-t Y s^/4 z i’/l TWpt/f (26) 


(with B — ^ - 1 ) and 


?• , tr-(z- zj 1 ] 




(27) 


while 


H= l 

r 0 


U X, -X >(( Y-V.it (Z-2J] (28) 

^ ^ - x £ [(X- Y)%(z-Z /]' /2 


and 


with 


[ r -= i 

r i 

L J f 

L 




T ( -- T . ©1 

M (*< - X ) ± R' 


(29) 


(30) 


* If the trailing edge is not fully supersonic then 
the contribution of the wake can be treated similarly 
to the subsonic case, with the only difference that the 
device of truncating the wake at finite distance is not 
necessary in the supersonic case, since only a finite 
portion of the w T ake can have an effect on the aircraft. 
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Using Eq. (12) and following the same procedure used 
for the subsonic case one obtains the supersonic indicial 
aerodynamic operator 

2 E(pJ T) -- Z B l ' ' (T- )+f( r ' ©')]+■ 

i c;[ 4 u (t-©:) + (31) 

E b \A- 4 > (t-© + ) + 4 , 4 > (T- ©;}] 

L 1 d l L ' J 

where 

3 ' ^ - -L $> M (p) JL dZ B 
c; , ©L | NjriZflU Zb 

^7i ^2 a L z?fj c i S'-' 

(32) 

!>' ^ _ J_ $> ft (?) JL ill dZ s 

t ^ *■ F In b 


In particular if P* coincides with the node J, Eq. (31) 


reduces to 


tf r K ^ ) * t' o ®« )] 

I K. i (T -^] 


where 


** - 

Finally, taking the Laplace transform of Eq. (33) one 


obtains 


m 
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where 



Ecruation (35) indicates that the matrix A ! is the 

t lui 

desired supersonic matrix transfer function relating 
the transformed vector of the ootential f i ! to the 

L A j ) 

transformed vector of the conormal-wash 1 Xl/' ? , for 

L *1 

supersonic trailing edge configurations. If the trailing 
■edge is not fully supersonic, the formulation may be 
modified following the same ideas, used for the subsonic 
case. It may be worth noting that according to Eq. (29) 
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CONCLUDING REMARKS 


A general theory for unsteady compressible potential 
aerodynamics has been presented. The motion is assumed 
to consist of small perturbation starting at time t=0) 
around a steady-state constant-velocity motion. In this 
case the relationship between the velocity potential 
and the normal-wash is given by an integral operator in 
space and a differential-difference operator in time. 

Using the finite element method to solve for the spacial 
problem one is left with a differential-difference equation 
in time. This can be solved numerically for time-domain 
or by using the Laplace transform and thus obtaining the 
matrix transfer function for frequency-domain analysis. 

The results presented here represent a considerable 
improvement with respect to the formulation available 
thus far since complex configurations could be analyzed 
only for steady and oscillatory flows (see, for instance, 
Refs, 3 and 4) while unsteady flows could be analyzed only 
for simple configurations such as zero-thickness wing with 
special planforms. With the method presented here unsteady 
flows around complex configurations can be analyzed for 
most cases of practical interest. For, the linear equations 
of flight dynamics implies small perturbations around a 
steady state motion, usually constant speed horizontal 
flight. The unsteady aerodynamic analysis presented here 
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does not require any additional limitation and therefore 
is the most general formulation within the above framework. 

In addition it should be noted that the increase 
in generality of the formulation is obtained at no 
additional increase in computational complexity. For, 
if simply oscillatory problems are considered the only 
advantage is to replace ^ with ifl (Ref. 3 and 4), 
with no particular computational saving. 

Another advantage of the present approach is that 

the dependence of the. matrix A JM upon the complex frequency, 

s, is given in a very simple explicity analytic form. 

This is a considerable computational advantage since once 

the coefficients B T . A , , D_ T , F_ M , G TVTf S„ r , and have 

LM JL JL JN JN NL JN 

been evaluated (and these are necessary even for the 
evaluation of the potential at one single reduced frequency) , 
it is a trivial matter (essentially the inversion of one 
matrix) to obtain the results at different values of s . 

Also advantage can be taken of the analytic dependence 
upon s to obtain approximate expressions for the matrix 
I • For instance, at low frequencies a Taylor series 

expansion for Eq. (24) can be taken to yield 
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where r -*-1 

[r,J=[^l-[ c .«!-[ F J[ s -l os,, 

L J 

It is essential to note the difference between the 
method presented here and the classical approach for 
unsteady (oscillatory) aerodynamics. In this case the 

solution is assumed to be oscillatory (i.e. of the type 

iaT , 4 

e i and then the problem is solved in space. Here the 

problem is solved first in space and then in time. This 
inversion of the time and space solutions might appear to 
be irrelevant, but is not. For, in the classical formula- 
tion the convergence of the space solution is analyzed on 
the time-transformed unknown which is highly oscillating 
in space. On the other hand, here the finite-element 
method is applied to the untransformed equation, where 
the unknown is smooth and therefore 
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fewer elements are required for convergence. The transform 
applies to the discrete system and therefore high frequency 
components do not involve a change in the number of elements. 
This question is analyzed more in detail in Appendix A, 
where the Laplace transform is used first and then the 
transformed equation is solved by finite elements. This 
process yields Eq. (A. 3} which should be compared to Eq, 

C20). Note that the integrals in Eq. (A. 4) are not frequency 
independent and therefore a large number of elements is 
required for convergence. On the other hand the integrals 
in Eq. (18) are frequency independent and the convergence 
must be evaluated on the time-domain solution. 

Another advantage of the present approach is the 

already mentioned possibility of truncating the wake at 

finite distance if the analysis is limited to T £ T max * 

Csince A&sO for T < T v <[],)• This eliminates the 

max M 

problem of convergence as the length of the wake goes 
to infinity. 

Other questions which have not been discussed here 
are the Kutta condition and the role of the diaphragms 
in supersonic flow. These points are analyzed in Refs. 

3 and 4 for a zeroth-order finite-element solution (i.e. 
potential, 4* / and normal-wash, ^ , constant within each 

element) . Further investigations for higher-order 
solutions are now under wav. 
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In conclusion a new approach for small-perturbation 
unsteady aerodynamics has been presented* Taking full 
advantage of the finite element method and the operational 
calculus the problem is simplified considerably and the 
relationship potential-normalwash is reduced to a system 
of algebraic equation, with explicit dependence of the 
coefficients upon the complex frequency, s . 
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APPENDIX A 


In this Appendix it is shown how interchanging 
the order of the time and space solution, a different 
type of equation is obtained. Taking the Laplace 
transform of Eq. (9) one obtains 

5 © 


4n £\'P,:g(Pj - - $ i? 


R 




Next using the finite -element method, i.e. setting 

L- 1 — 

#fPh 2 4 Nj?i 

aJ(P) r 2 L^P) 


A. 2 
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yields, at P - P L , 


= [s IL ] { + L ^ J Hu \ + 5 fA*. ] Hu i 

PjN I i> AC h.j ^ 4 S L^JW i i ^ j 


A, 3 


where 


B. 


JL 


_ _L £> M ( p ) e 

. ** "i. 


- S 


e 


R 


6^i ^ p 

G> 


Z * ?J 


C = i JL £t Njp) e ' S 3 — (-S-) _ 

Jl Unit u ' ZtJ'Ri J?,= P, 


|/*> 


B 


t>. 


T " 


3L 




.1 <S Njp) «-‘ S i §2 -< 2„1 

-2/1 ^Z B R -iP,= ? J 

— ff L (P) e %& A. fJ. )d2 w 
2 ” JJ iw N <?aAR' w 


G 




J_ ff 


, fp) *- s0 _L 25 d2, 
l m' k, - € r aro w 


p* u p 5 

A. 4 


P s 1> 

* j 


J;7 JJ y 

d. ^ 

Finally, according to Eqs. (15) and (16) one obtains 

- - S 77 ^ r . - s 7 }, 


/ ^ ^ _ s // r 



Combining Eqs. (A. 3) and (A. 5) ahdf solving for £ <rp^ £ 

yields a new expression for the matrix transfer function 
r 




\f +s 5 Ifs e ' sn "1 
i 'jrt 3N J L j 


j \* e 


- $ S t 


lh 


A. 6 
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Comparison with Eq. (24) indicates clearly the advantage 
of using the space-discretization before the time-trans- 
formation. 
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APPENDIX B 


In this Appendix Eq. (10) is derived. Consider 
the Bernoulli theorem for potential flow 


M , i ifS/ z + f is. -- 

dt 2 1 ' * J f 


CB.l) 


Since no pressure difference can exist across the wake, then 


(fowl'd F*i !v <> =° 


(B. 2 } 


~ ‘fc) + i pt- 7 k)‘W t ■ H)y ° 


This can be rewritten as 




CB.3) 


(B. 4 } 


jD t 


(%- <h)‘° 


where 


= i. + 


(B.5) 


(B. 6) 


is the total time derivative obtained by following 
a particle having the average velocity. 


J (f</> * 


CB.7) 
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Equation (B .5) implies that the wake vortices are 
simply convected, that is that 

' * f ' ;7 ) (B - 8) 

where XI is the time necessary for a point of the wake 

to move from a point P TE of the trailing edge to the 

point P. If the unsteady flow is infinitesimal 71 can 

be evaluated from the steady state solution. If small 

perturbation apply to steady state flow as well, n is 

given by the distance i: ~* TE divided by the velocity. 


51 = ^ 


(B , 9) 


Using nondimen sional variables (Eq. 2), Eqs. (B.8) and 
CB.9) are replaced by Eqs. (10) and (11) respectively. 
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